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Abstract 

Bogomol'nyi-Prasad-Sommerfield (BPS) vortices in U{N) gauge theories have two layers 
corresponding to non-Abehan and Abehan fluxes, whose widths depend nontrivially on the 
ratio of U{1) and SU{N) gauge couplings. We find numerically and analytically that the 
widths differ significantly from the Compton lengths of lightest massive particles with the 
appropriate quantum number. 
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1. Introduction. Many important properties of Abelian (ANO) vortex were found [H El El SI E] 
since its discovery [6]. Recently vortices in U{N) gauge theories (called non- Abelian vortices) 
were found [71 [8] and have attracted much attention ^ because they play an important role in 
a dual picture of quark confinement [HI [ID] and are a candidate of cosmic strings [11] (see [12] 
for review report). The moduli space of U{N) non- Abelian vortices was determined in [13] and 
study on interactions between non-BPS configurations started in [TJ]. Non- Abelian vortices in 
other gauge groups have been studied 



Although there have been much progress and wide applications, internal structures and de- 
pendence on gauge coupling constants have not yet been studied for (color) magnetic flux tubes. 
It is particularly important to study physical widths of vortices qualitatively and quantitatively, 
although it is not easy because no analytic solutions are known. It may be tempting to speculate 
that the width is determined by the Compton lengths of lightest massive particles with the ap- 
propriate quantum number. Purpose of this letter is to clarify intricate multiple layer structures 
of non-Abelian vortices by investigating numerically and analytically the equations of motion. 
Non-Abelian vortices have two distinct widths for SU{N) and U{1) fluxes. We clarify properties 
of these widths by making use of several approximations. It turns out that non-Abelian vortices 
are very different from ANO vortices and have much richer internal structures. 

2. Vortex equations and solutions. Let us consider a U{N) gauge theory with gauge fields 
for SU{N) and for U{1) and Higgs fields H {N-hj-N matrix) in the fundamental repre- 
sentation. We consider the Lagrangian C = K — V which can be embedded into supersymmetric 
theory with eight supercharges 
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K = Tt 



2g 



4e2 



q2 . .2 
V = — Tr 



^__^{HH^y]+'-{TT[HH^-cl^]Y, (2) 

where {X) stands for a traceless part of a square matrix X. Our notation is T>nH = {d^ + iW^ + 
iw^lN)H, F^^ = d^W^ - dyW^ + i \W^, W^] and f^^ = d^w^ - d„w^. We have three couplings: 
SU{N) gauge coupling g, U{1) gauge coupling e and Fayet-Ihopoulos parameter c > 0. 

The Higgs vacuum H = ^/clN is unique and is in a color-flavor SU{N)c+f locking phase. 
Mass spectrum is classified according to representations of SU{N)c+f as = gy/c for non- 
Abelian fields 0N = (W^5 {H)) and me = eV2Nc for Abehan fields 0a = (ui, Tt{H — y/clN)). The 
non-Hermitian part of H is eaten by the U{N) gauge fields. A special case of mg = me [7] has 
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been mostly considered so far, which is equivalently 

7 = 1 with 7 = ^= = 3 

but we study general cases in this Letter. 

Let us consider static vortex-string solutions along xs-axis. The BPS equations for the non- 
Abelian vortex are 

2c rftg 2c 

with P = ("Di + iD2)/2. The tension of A;-vortex is = — c j d?xTi[fi2lN] = ^irkc. No 
analytic solutions have been known whereas a numerical solution was first found in [8]. For 
k = 1 vortex, we take H = S-^Hq, W = -iS'^dS [2W = {Wi + wiIn) + ^(Ws + wsl^v)] with 
diagonal matrices Hq = diag( re^^ 1, ■ ■ ■ , 1 ) and cSS^ = ei^^+^^°^''')^^+i^^+^^''^''')^ . Here, 
T = diag(l, — ]y3y, ■ ■ ■ , ~7v^) and ipg are real functions of the radius r > of the polar 

coordinates (r, 6*) in xi,X2 plane. Then we get 



■>l>q \ 

+ (A'-l)e— j =1, (5) 
^"i + ^e-'-fe-^-eA^-O, (6) 



Aibe 1 

mi JM 



ml N 

with A/(r) = dr{rdrf{r))/r. The boundary conditions are ■j/'e, V'g ('^ oo) and N^e, j^Zii^g ~ 
— logr^ (r ^ 0). The fluxes and the Higgs fields are expressed by 

/i2 = -^AV'e, Fi2 = -^AV'gT, if = ycdiag(/ii,/i2,- ■■ ,^2) (7) 

with hi = e 2 i"**^ and /i2 = e~ 2 ]v-i>. xhe amount of the Abelian flux is and the 
non-Abelian flux is (A^ — 1)/A^ of the ANO vortex. 

We found the following theorems for Eqs. ([5]) and 

a) 'ipcg > 0, dr4'e,g < and A4'e,g > 

b) < |/i2|, I /ill < 1 and 9^1 > 

c) dr\h2\ I and 1 I /i2 I ^/N/{P^^f^ for 7 | 1 

All these can be proved by using the following theorem for an analytic function /(r) satisfying 
/(r) < ^ A/(r) < 0: If5J(0) < Oand/(cx)) = 0, then /(r) > OforV G (0,oo). In the case of 
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Fig. 1: hi^2 (solid, broken lines) and i?e,g- The left panels [rrie = 1) for log7 = 0, 1, 2, 3, 4, 5 and oo. 
The right panels {nig = 1) for log7 = 0, —0.5, —1, ■ ■ • , — 3 and — oo. 

7 = 1, we get Nipe = j^zii^g = "^ano and the above equations reduce to At/^ano = 'ml{l — e~'^^'^°) 
with boundary condition ipANO ^ (r — » oo) and ipANO ~^ — logr^ (r — > 0). 

Numerical solutions for N = 2 for a wide range of 7 (including 7 = 0, 00) are shown in 
Fig. [TJ Winding field hi is not sensitive on 7 while unwound field h2 is. As rrig being sent to 00 
(7 —>■ 00), the non-Abelian flux F12 becomes very sharp and finally gets to singular. Interestingly, 
the Abelian flux /12 is kept finite there. In a region 7 < 1 (me > rrig), on the other hand, the 
Abelian flux is a bit smaller than the non-Abelian tube. Surprisingly, the fluxes remain finite 
even in me 00 limit. 

3. Asymptotic width. Let us investigate the vortex solution by expanding ([5]) and in region 
r ^ max{m~ m"-*^} where |'?/'e|; I'^gl ^ 1- We keep only the lowest-order term in ip^ while 
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keeping terms up to next to leading order in ipg in Eq.(l5]): 

(A - ml)^e + ^^/^ = 0, (A - mDi^g = 0. 
The solution is given by the second modified Bessel function Ko(r), and approximated as 



'■\/ 2me 



9 — ^9 



TT 



2mgr 



(9) 



4(Af-l){l-472) nigr ' 

with Ce,g being dimensionless constants, see Table [H The asymptotic behavior of ipe changes at 
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Ce 


C9 


be 


b9 


a. 







1.1363 





0.75905 


V2 


0.25 




1.1853 


0.31719 


0.73163 


1.31688 


0.5 




1.3090 


0.47907 


0.68393 


1.18361 


0.75 


2.196 


1.4852 


0.55921 


0.64006 


1.07932 


1 


1.70786 


1.70786 


0.60329 


0.60329 


1 


1.5 


1.4715 


2.3031 


0.64726 


0.54697 


0.88820 


2 


1.4037 


3.15 


0.66773 


0.50604 


0.81226 


2.5 


1.3746 


4.32 


0.67897 


0.47469 


0.75640 


3 


1.3594 


6.0 


0.68584 


0.44969 


0.71301 


00 


1.3267 




0.70653 









Table 1: Numerical data for = 1 f/(2) vortex. 

7 = 1/2 (upper for 7 > 1/2 and lower 7 < 1/2). Similar phenomenon was observed for the non- 
BPS ANO vortex (3^|4j. The origin of %l)e{ipg) is (non-)Abelian fields 0a(0n) with mass meimg), 
and the 7 = 1/2 threshold can be interpreted as follows. The expansion of the Lagrangian with 
respect to small 0a,n contains the triple couplings 0a0n- For nie < 2mg, asymptotics for 0a,n 
are given by KQ{mf.^gr) as the two-dimensional Green's function. When mg > 2mg, the particles 
decay into two particles 0^ through these couplings, and thus, (px exhibits the asymptotic 
behavior e"^™^'' below 7 = 1/2 like Eq. ([9]). On the contrary, even for 7 > 2, 0n does not behave 
as e~^™"='' since there is no triple coupling 0n0a "^^^ traceless condition for 0n- 

Let us define asymptotic width of the vortex by an inverse of the decay constant in Eq.Q: 

2 /me for 7 > 1/2, 

L 

2/{2mg) for 7 < 1/2, 
4 



2/mg 



(10) 



Here the factor 2 is put in the numerator to match with another definition in Eq.( |T4l) . The 
asymptotic width of Abehan vortex is bigger than the non-Abehan one when 7 > 1 and vice 
versa for 1/2 < 7 < 1. For 7 = 1, the two widths are the same. The case 7 < 1/2 indicates a 
significant modification, where the Abehan fiux tube is supported by the non-Abehan fiux tube. 
This answers the question why the Abehan vortex does not collapse even in the 171^ —>■ 00 hmit. 
When 7^1, the thin non-Abehan fiux hidden by fat Abehan fiux cannot be correctly measured 
by Lg. We now turn to another definition of vortex width which refiect the size of the vortex 
core more faithfully. 

4. Core widths. Let us consider a region near the vortex core. We expand fields by 



tpe^-^^ogKrrier, i)g^-2^^^-j^\oghgmgr. (11) 



Dimensionless constants fee,g are related to h2{r = 0) by 

a, = hiO) = ihehbgY"". (12) 

See Table [H These are important since they are related to the maximum values of the magnetic 
fiuxes at r = 

5e--— l^l-^^Sj, i?.--^^^«r (13) 

Widths of the magnetic fiuxes can be estimated by using a step function Q{x) as F12 = BgQ{Lg — 
r)T and /12 = Bf.Q{L^ — r) keeping amounts of the fiuxes as |-Be| x vrLg = 2tt/N and x 71 = 
2{N-l)n/N: 

2 ~ 2 

Le = , Lg = , (14) 



We call Le and Lg as the core widths of the vortex. In the case of 7 = 1, Le and Lg coincide 
because of a^=i = /12 = 1. In Fig. [21 we show the core widths numerically in the case of = 2, 
which are analytically reinforced as we will discuss. We again observe that the Abelian core does 
not collapse even when mg ^1 (7 <^ 1). 

Mass dependence of the core widths coincides with one of the asymptotic widths L^^g given in 
Eq. ffTOj) . except for I/g(7 > 1), see Fig. [2l The asymptotic width Lg is independent of me whereas 
the core width Lg depends on nie- This is because Le,g more faithfully refiects the multilayer 
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structure in the large intermediate region of r for the strong coupHng regime (7^1), to which 
we now turn. 

5. Two strong coupling limits. Here we study two hmits: i) — > 00 with mg fixed, and ii) 
rrie 00 with irig fixed. In the former hmit (7 — > 00 with me fixed), all the fields with mass 
mg become infinitely heavy and integrated out from the theory. As a result, the original U{N) 
gauge theory becomes the Abelian theory with one SU (N) singlet field B = det H. Note that 
the target space is C/Zn because U{1) charge of B is A^. Eq. ([6]) is solved by '4'g,oo = while 
'ijje,oo is determined by Eq. ([5]) 

= 1 - e-'^-"", Ve.oo -^log^e&e.oor (15) 

mj r-*o N 

where suffix oo denotes 7 = oo: ipg^oo = i'gl-y^oo- The boundary condition tells that vorticity is 
fractional k = 1/N. This way the non- Abelian flux tube collapses and the U{N) non- Abelian 
vortex reduces to the 1/A^ fractional Abelian vortex. This solution helps us to understand the 
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non-Abelian vortex for 7^1 better. Since e ^ (^e&e.oo''')^^^ <^ 1/me) for j 1, ipg for 
r -C 1 / TUf, is well approximated by a solution of the following 

^ + ^Kr)^ fe-^« - e^) = 0, (16) 
rrig N \ / 

where the parameter rhg = 'n^g(&e,oo7~^) has a mass dimension. Therefore ipg has asymptotic 
behavior in the middle region 1 / rhg <^ r ^ 1/ rrie 

and ^/'g ~ — ^^^^^ log (6g m^r) for r ^ ^/'r^g- Here 6^, are determined numerically and inde- 
pendent of 7, for instance, bg = 0.74672, Cg = 0.63662 for N = 2. Comparing this with Eq. f|TT]) . 

1 1 
we find bg ^ bg [6e,oo7~^]'^ and ^bg^ \be,ool~^\^^ for 7 > 1- 



In the second limit (7-^0 with nig fixed), all the fields with the mass m-e are integrated out. 
The model reduces to a CP^^~^ model with SU{N) isometry [in SUlN"^ — 1)] gauged. Eq. 
is solved by 6^"'° = (e~'^9'° + (A^ — l)e'^-'^)/N while ip9,o is determined by 

{N - 1) (1 - e^^^^'A 
^^9,0 = ml— ^ ^— ^, (18) 

where the suffix denotes 7 — > 0: ipg^o = ipgl^y^Q. This is a new a-model lump with the non- 
Abelian fiux accompanied with the internal orientation CP^^^. Again we can make use of this 
solution to understand the non-Abelian vortex for 7 ^ 1. Let us define = "^^^'"1°^ = ^1 
which turns out to be finite a'^ = {N ~ 1)/(1 + 46^|'y_,o)- Since Aipg Q^O) = and Aipe,o{0) = 
\im^_,oml-f-\l - {N - l)a^^/N), we find = ao (1 - iVC^a^) + ■■■)> «o = V^/i^ " 1) ^r 
7 < 1. 

6. Summary and Discussion. We have proposed two length scales for fiuxes of non-Abelian 
vortices: asymptotic widths L^^g in Eq. ffTOj) and core widths L^^g in Eq. f|T^ . By using the 
asymptotics of obtained above, the core width is summarized as 



2a 2 /Af-1 



\/¥} (7 « 1), 



where a and (3 depend only on A^ and are determined numerically, for instance a = 0.55010, j3 = 
0.97022 for N = 2. The core and asymptotic widths have the same mass dependence except 



for Lg and Lg for 7^1. Interestingly, the Abelian flux does not collapse even when rrie 00 
(7 <^ 1). For 7^1, the thin non- Abelian flux is hidden by fat Abelian flux, so that the true 
width cannot be captured by the asymptotics at r ^ ^e^- Instead we should use improved 
approximation given in Eq. (ITTl) to measure the non- Abelian flux. Indeed, the decay constant in 
Eq. f|T71) is fh~^ whose mass dependence is the same as one of Lg for 7^1. 

In the limit — > 0, the original U{N) gauge theory reduces to U{1) gauge theory coupled 
to A^^ Higgs fields. Eq. f|T9l) tells us that the vortex is diluted and vanishes in this limit. This 
is consistent with the fact that there is no (smooth) vortex solution with a winding number 
1/A^ in that U{1) theory. The minimal vortex in the U{1) theory corresponds to A^ vortices in 
the original theory. The dilution is expected to be avoided and all the fields with mass of the 
order of irig decouple, when A^ vortices are arranged as H = f{r)l]\[. In the limit mg — > 0, 
there is no BPS vortex solution since the U{1) gauge field is decoupled from the Higgs fields. 
Actually, according to Eq. f|T9|) . one can find both of the Abelian and the non- Abelian fluxes 
are diluted again even in this limit due to the factor 7^+1 . Monopoles/instantons attached by 
vortices are known to exist [10]. The above observation implies that such configurations reduce 
to a monopole/instanton configurations in the SU{N) gauge theory in that limit, and strongly 
supports the correspondence between the moduli spaces of them. 

It is interesting to study relation between non-BPS vortices and monopoles. It was found 
that monopoles do not collapse when the Higgs mass is very large [16]. 
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